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1  INTRODUCTION 


Millimeter  wave  dielectric  waveguides  have  been  extensively  studied  during  the  past  two  decades. 
Examples  of  these  waveguides  include  dielectric  image  guides  [1].  strip  dielectric  guides,  insulated 
image  guides,  strip-slab  guides  [2],  inverted  strip  dielectric  guides  [3].  cladded  image  guides  [4]  and 
trapped  image  guides  [5].  These  waveguides  are  constructed  from  combinations  of  layers  and  ridges 
of  various  permittivities  in  order  to  provide  a  region  wherein  the  propagating  power  is  well-confined. 
The  widths  of  these  lines  approach  one  guided  wavelength  in  order  to  maximize  field  confinement. 
Although  there  are  several  examples  of  the  monolithic  use  of  these  dielectric  lines  in  the  literature 
[6],  they  have  been  generally  considered  hybrid  in  nature. 

Recently,  Engel  and  Katehi  [7]  suggested  the  development  of  monolithic  sub-mm  guiding  struc¬ 
tures  that  can  be  realized  by  considering  variations  of  the  early  dielectric  lines.  The  new  waveguides 
are  made  of  materials  which  are  available  in  monolithic  technology  so  that  they  are  compatible  with 
solid-state  sources  (see  Fig.  1).  The  dimensions  of  these  monolithic  guides  are  fractions  of  a  guided 
wavelength,  so  the  new  structures  may  be  used  not  only  as  guiding  media  but  as  means  of  making 
passive  components. 

The  successful  use  of  dielectric  lines  in  a  hybrid  or  monolithic  environment  relies  mainly  on 
the  ability  to  realize  an  efficient  transition  to  the  dielectric  waveguide.  In  the  past,  a  variety  of 
transitions  from  rectangular  waveguide  to  mm-wave  dielectric  lines  have  been  mainly  characterized 
experimentally  [8,  9].  When  layered  ridge  dielectric  waveguides  (LRDWs)  [7]  are  used  in  a  monolithic 
environment,  transitions  to  and  out  of  dielectric  waveguide  may  be  realized  through  a  short  length 
of  printed  conductor  line  on  top  of  the  ridge  (Fig.  2)  [10,  11].  In  this  manner,  the  dielectric 
waveguide  components  are  effectively  coupled  to  other  conventional  monolithic  circuit  components 
which  are  printed  on  the  same  wafer.  Previous  work  on  the  characterization  of  such  a  transition 
exclusively  considered  the  transition  operating  in  an  ideal  shielding  environment  which  was  designed 
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Figure  1:  Cross-section  of  the  layered  ridge  dielectric  waveguide  (LRDW)  structure  under  study. 

to  effectively  shield  only  that  circuit  component.  As  a  result,  the  cavity  was  chosen  adequately 
small  enough  and  it  was  placed  far  enough  from  the  transition  so  that  it  did  not  resonate  within 
the  frequency  band  and  did  not  interfere  with  the  circuit.  The  electrical  performance  of  this  ideal 
transition  was  analyzed  by  a  hybrid  full-wave  integral  equation-mode  matching  (lEMM)  analysis 
technique  and  preliminary  results  were  presented  in  [10].  Furthermore,  a  detailed  study  of  the 
same  shielded  transition,  using  both  the  lEMM  and  the  Finite  Difference  Time  Domain  (FDTD) 
methods,  is  presented  in  [11].  In  practice  however,  monolithic  circuits  operate  in  open  environments 
or  within  a  larger  shielding  packages  and  may  suffer  from  parasitic  radiation  (space  and  surface  waves 
losses)  which  limits  performance  considerably  and  intensifies  unwanted  electromagnetic  interference 
between  circuit  components.  As  a  result,  parasitic  radiation  needs  to  be  studied  carefully  and  should 
be  accounted  for  during  design.  In  this  report,  the  above  transition  is  studied  in  open  environment 
using  the  FDTD  technique  and  the  effects  of  parasitic  radiation  on  the  transition  performance  are 
discussed.  Moreover,  the  transition  from  a  conventional  mirostrip,  printed  on  the  substrate,  to  the 
LRDW  is  analyzed.  The  possibility  of  using  other  transitions  is  also  discussed. 
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Figure  2:  Strip-ridge  line  to  LRDW  transition.  Strip  \vidth=:20  /im 

The  second  part  of  this  work  is  devoted  to  the  dielectric  waveguide  coupler.  The  distributed 
directional  coupler,  shown  in  Fig.  3,  which  consists  of  a  section  of  a  coupled  dielectric  guide  of 
length  L  and  separation  S,  presents  the  simplest  coupler  geometry.  Such  a  coupler  can  be  made 
of  any  type  of  dielectric  waveguide  and  has  been  the  subject  of  study  of  many  publications  (see 
for  example  [12]-[14]).  In  most  of  these  studies  the  coupler  design  was  completed  by  choosing  the 
appropriate  coupling  length  through  empirical  formulas  for  the  magnitudes  of  Si2  and  S13.  This 
approach,  however,  neglects  interference  between  the  feeding  lines  as  well  as  radiation  loss  which  can 
be  substantial  in  open  dielectric  waveguides  [15]  and  predicts  performances  which  are  unrealizable. 
To  overcome  this  problem,  when  comparing  experimental  and  theoretical  results,  measured  data 
had  to  be  normalized  accordingly  to  numerically  eliminate  the  effect  of  the  dielectric  and  radiation 
losses  at  the  bends  and  junctions  [13,  14].  In  this  report,  the  FDTD  technique  is  used  to  analyze  the 
coupler  geometry  and  simulates  the  open  environment  in  which  the  coupler  operates  by  appropriate 
absorbing  boundary  conditions.  For  validation  purposes,  a  mm-wave  coupler  which  utilizes  the 
image  guide  is  analyzed,  first,  and  compared  to  experimental  data.  Following  this  validation,  a 
sub-mm  wave  directional  coupler  made  of  the  LRDW  shown  in  Fig.  1  is  extensively  studied.  The 
scattering  parameters  and  radiation  loss  factor  of  both  couplers  are  presented. 

The  third  part  of  this  work  concerns  the  issue  of  matching  the  LRDW.  Looking  at  the  field 
distribution  on  a  cross  section  of  the  LRDW  indicated  that  most  of  the  power  is  confined  in  the 
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Figure  3:  Top  view  of  a  distributed  directional  coupler. 

substrate  just  underneath  the  ridge.  The  ridge  acts  only  as  a  guiding  element  for  the  power  propa¬ 
gating  along  the  line.  With  this  fact  in  mind,  an  open-end  LRDW  (by  stopping  the  ridge  at  some 
point  while  leaving  the  substrate  extended)  is  analyzed  and  the  return  loss  is  computed  using  the 
FDTD  technique.  It  is  shown  that  such  an  open-end  may  be  used  as  a  matching  termination. 
Alternative  matching  terminations  are  also  discussed. 

2  THEORY 

The  FDTD  method  was  first  introduced  by  Yee  [16]  to  solve  electromagnetic  scattering  problems. 
In  this  method,  Maxwell’s  curl  equations  are  expressed  in  discretized  space  and  time  domains  and 
are  then  used  to  simulate  the  propagation  of  an  initial  excitation  in  a  “leapfrog”  manner.  Recently, 
the  method  has  been  successfully  applied  to  characterize  microstrip  and  coplanar  waveguide  (CPW) 
lines  and  discontinuities  [17]-[24],  optical  integrated  circuits  [25,  26],  and  mm-wave  dielectric  line 
transitions  [27].  Only  a  brief  summary  of  the  FDTD  method  is  described  here. 

In  order  to  characterize  any  planar  discontinuity,  propagation  of  a  specific  time-dependent  func¬ 
tion  through  the  structure  is  simulated  using  the  FDTD  technique.  The  time  dependence  of  the 
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excitation  can  be  chosen  arbitrarily.  A  Gaussian  ])ulse  is  olten  used  because  it  is  smoothly  varying 
in  time  and  its  Fourier  transform  is  also  a  Gaussian  function  centered  at  zero  frequency.  However, 
since  dielectric  lines  possess  a  non-zero  cutoff  frequency,  it  is  more  convenient  to  use  an  excitation 
pulse  which  has  an  amplitude  spectrum  in  a  specific  frequency  range.  In  this  research,  the  following 
function  is  used  [28]: 


m 

T 


(‘-‘cr  _ 

Ce  ~2t2  sinfief) 

n/6 

^{fmax  fmin) 


(1) 

(2) 


^  —  '^{fmax  d“  fmin) 


(3) 


to  —  3T 


(4) 


where  fmin  and  fmax  are  the  frequencies  at  which  the  value  of  the  spectrum  is  5%  of  its  peak 
which  occurs  at  {fmax  +  /mm)/2.  The  above  function  is  the  imaginary  part  of  what  is  known  as 
the  Gabor  elementary  function  [29].  Following  the  time  and  space  discretizations  of  the  electric 
and  magnetic  field  components,  the  FDTD  equivalents  of  Maxwell’s  equations  are  then  used  to 
update  the  spatial  distributions  of  these  components  at  alternating  half  time  steps  [30].  The  space 
steps,  Ax,  Ay  and  Az,  are  carefully  chosen  such  that  integral  numbers  of  them  can  approximate 
the  various  dimensions  of  the  structure.  As  a  rule  of  thumb  and  in  order  to  reduce  the  truncation 
and  grid  dispersion  errors,  the  maximum  step  size  is  chosen  to  be  less  than  1/20  of  the  smallest 
wavelength  existing  in  the  computational  domain  (i.e.,  at  the  highest  frequency  represented  in  the 
pulse).  Then,  the  Courant  stability  criterion  is  used  to  select  the  time  step  to  insure  numerical 
stability.  For  the  LRDW  structures  under  study,  the  following  parameters  are  used:  Ax=5  fim, 
Ay=Az=lO  fim  and  At=0.0136  ps. 

In  order  to  excite  the  dielectric  waveguide,  the  vertical  electric  field  component  at  the  front  plane 
(z=0)  is  excited  (see  Fig.  4)  and  the  magnetic  wall  source  condition  of  [18]  is  used  to  compute  the 
fields  elsewhere  in  the  plane  z=0.  The  source  distribution  has  been  modified  to  take  into  account 
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Figure  4:  Excitation  for  the  dielectric  waveguide  in  the  FDTD  analysis. 

the  discontinuity  experienced  by  the  vertical  electric  field  [19].  For  the  electric  field  components 
lying  on  a  dielectric-dielectric  interface,  the  average  between  the  two  permittivities  is  used  in  the 
FDTD  equations  [17].  The  super-absorbing  first-order  Mur  boundary  condition  [31,  32,  33]  is 
utilized  to  terminate  the  FDTD  lattice  at  the  front  and  back  planes  in  order  to  simulate  infinite 
lines.  This  absorbing  boundary  condition  requires  a  choice  for  the  incident  velocity  of  the  waves, 
or  equivalently  Ce//-  It  has  been  found  that  an  appropriate  choice  of  Cr.e//  minimizes  the  effect 
of  the  absorbing  boundary  walls.  For  the  dielectric  waveguide  under  consideration,  the  average 
between  the  relative  effective  dielectric  constants  at  and  fmax  is  used.  These  effective  dielectric 
constants  can  be  obtained  from  the  2D-FDTD  method  [34,  35]  or  by  simulating  a  propagating  pulse 
on  a  dielectric  through  line  [36].  The  effectiveness  of  these  choices  has  been  checked  by  performing 
several  numerical  experiments  with  different  values  of  Cr.e//-  Oii  the  other  hand,  the  first-order 
Mur  boundary  condition  is  used  on  the  top  and  side  walls  to  simulate  an  open  structure.  For  the 
dielectric  waveguide  under  study,  the  left  and  right  absorbers  were  placed  at  a  distance  of  240  /um 
from  the  edges  of  the  dielectric  line,  while  the  top  absorber  was  placed  at  a  distance  of  300  //m  from 
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the  top  surface  of  the  line. 


In  general,  the  frequency  dependent  scattering  parameters,  8,^,  can  be  obtained  as  follows  [17,  18]: 

,, ,  ,  y.i-  ~ 


I  Zqj 
Zoi 


(5) 


where  V,  and  Vj  are  the  voltages  at  ports  i  and  j,  respecti\'el\’.  and  Zqj  and  Zoi  are  the  wave 
impedances  of  the  lines  connected  to  these  ports.  In  case  of  dielectric  lines,  where  voltage  is  not 
uniquely  defined,  the  vertical  component  of  the  electric  field  is  probed  and  used  in  the  above 
equation  instead  of  the  voltage  [15].  To  obtain  Sn(a;),  the  incident  and  reflected  fields  must  be 
known.  Since  the  FDTD  simulation  calculates  the  total  field  which  is  the  sum  of  the  incident  and 
reflected  waveforms,  the  incident  field  is  obtained  from  a  line  of  an  infinite  extent  and  is  subtracted 
from  the  total  waveform  to  yield  the  reflected  field. 


3  RESULTS  and  DISCUSSION 

3.1  Transitions  to  the  LRDW 

In  this  section,  several  transition  structures  to  the  LRDW  are  studied.  Namely,  strip-ridge  to 
LRDW  transition  (Fig.  2),  conventional  microstrip  to  LRDW  transition  (Fig.  5)  and  rectangular 
waveguide  to  LRDW  transition  (Fig.  6)  [37]. 

Before  analyzing  these  structures,  one  needs  to  look  at  the  dispersion  diagrams  of  both  the  LRDW 
and  the  strip-ridge  line.  Figures  7  and  8  show  the  dispersion  curves  of  the  dominant  propagating 
modes  of  the  strip-ridge  line  and  the  LRDW,  respectively.  The  3D-FDTD  results  are  obtained  by 
simulating  the  propagation  of  a  pulse  through  an  infinite  line  [36].  On  the  other  hand,  the  2D-FDTD 
results  are  obtained  in  an  implicit  way  where  one  inputs  a  value  for  the  propagation  constant  and 
then  finds  the  frequencies  at  which  the  modes  propagating  with  this  specific  propagation  constant 
can  be  supported  [34,  35].  The  3D-FDTD  has  the  advantage  that  only  one  simulation  is  needed 'to 
obtain  the  dispersion  curve,  while,  the  2D-FDTD,  though  faster  and  simpler  to  implement,  needs 
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Figure  5:  Microstrip  to  LRDW  transition. 


9 


LRDW 


Figure  6:  Rectangular  waveguide  to  LRDW  transition. 
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Figure  7:  Dispersion  diagram  for  the  dominant  mode  of  the  strip-ridge  line  with  strip  width  of  20  fxm 
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Figure  8:  Dispersion  diagram  for  the  dominant  mode  of  the  layered  dielectric  line. 

to  be  executed  several  times  to  give  the  propagation  constant  in  a  frequency  range.  One  added 
advantage  of  the  2D-FDTD  technique,  however,  is  its  ability  to  predict  the  propagation  constants 
of  higher  order  modes  besides  the  dominant  propagating  mode.  Figure  8  shows  that  the  dielectric 
guide  dominant  mode  has  practically  a  cutoff  frequency  around  350  GHz. 

1.  Strip-Ridge  to  LRDW  Transition  (Fig.  2) 

Of  primary  interest  in  evaluating  the  performance  of  the  the  transition  shown  in  Fig.  2  is  the 
power  transferred  from  the  dominant  mode  in  the  strip-ridge  structure  to  the  dominant  mode 
in  the  LRDW.  Figure  9  shows  the  magnitude  of  the  reflection  coefficient  for  this  transition  in 
the  frequency  range  DC-550  GHz.  Due  to  the  fact  that  the  strip-ridge  to  LRDW  transition  is 
an  open  structure,  the  magnitude  of  is  always  less  than  1.0  even  for  frequencies  below  the 
waveguide  dominant  mode  cutoff  frequency  (around  350  GHz)  while  the  shielded  transition 
showed  a  total  reflection  (15111  =  1)  in  [11].  Specifically,  in  the  open  transition,  for  operating 
frequencies  between  100  and  300  GHz  almost  20%  of  the  incident  power  is  lost  in  the  form  of 
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Figure  9:  Magnitude  of  Sn  for  the  strip-ridge  to  LRDW  transition  shown  in  Fig.  2. 

space  and  surface  waves  which  will  be  referred  to  as  radiation  loss  in  this  report.  However,  as 
the  operating  frequency  approaches  500  GHz,  the  return  loss  reduces  to  around  -15  dB. 

2.  Microstrip  to  LRDW  Transition  (Fig.  5) 

Figure  10  shows  the  magnitude  of  the  reflection  coefficient  for  this  transition  in  the  frequency 
range  DC-550  GHz  for  different  strip  widths.  It  can  be  seen  that  the  return  loss  for  strips’ 
widths  of  40;um  and  60  pm  is  almost  the  same  and  slightly  better  than  that  obtained  with 
strip  width  of  20pm.  In  comparison  with  the  strip-ridge  to  LRDW  transition,  the  microstrip 
to  LRDW  transition  exhibits  a  larger  return  loss.  In  all  the  above  examples  the  dielectric 
waveguide  ridge  has  a  layered  consistency  similar  to  the  one  on  Fir.  1. 

3.  Rectangular  waveguide  to  LRDW  transition  (Fig.  6) 

Figure  6  shows  a  transition  from  a  recy‘tangular  waveguide  operating  in  the  dominant  mode 
from  350GHz  to  530GHz  to  a  LRDW  with  the  cross  section  shown  on  Fig.  1.  Figure  11  shows 
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Figure  10:  Magnitude  of  Sn  for  the  conventional  microstrip  to  LRDW  transition  shown  in  Fig.  5.  W  is  the  mirostrip 
width. 

the  magnitude  of  the  reflection  coefficient  for  this  transition  in  the  frequency  range  350-530 
GHz.  It  can  be  seen  that  the  reflection  coefficient  for  such  a  transition  is  less  than  14  dB  in 
the  whole  frequency  range. 

3.2  LRDW  Directional  Coupler 

A  simple  design  of  the  dielectric  directional  coupler  is  shown  in  Fig.  3,  where  9  indicates  the  angle 
between  the  tilted  line  and  the  horizontal  axis,  L  is  the  length  of  the  coupled  LRDW  sections,  s  is 
the  separation  between  the  coupled  sections  and  d  is  the  separation  between  the  feeding  lines.  This 
design  approach  assumes  zero  loss,  negligible  coupling  due  to  the  feeding  dielectric  waveguides  and 
the  tapered  sections  and  considers  quasi-static  interactions  in  the  coupling  region.  As  a  result,  the 
magnitudes  of  the  scattering  parameters  S12  and  S13  can  be  obtained  from  the  following  relations 
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Figure  11:  Magnitude  of  Sn  for  the  rectangular  waveguide  to  LRDW  transition  shown  in  Fig.  6.  The  waveguide 
has  a  width  and  height  of  560^m  and  280/im,  respectively. 


[12]-[14]: 

I  512  1  =  1  COS  ^ . —  l)j  I  (6) 

I  5i3  1=1  sin  2  I 

where  and  (3o  are  the  even  and  odd  propagation  constants  of  the  coupled  dielectric  waveguide, 
respectively.  With  the  knowledge  of  /?*  and  /3o,  the  coupling  length  L  can  be  chosen  to  provide  a 
coupler  with  a  specific  coupling  coefficient  jS'isI-  An  improvement  of  such  an  ideal  design  can  be 
achieved  by  extending  the  effective  coupling  region  within  the  tapered  sections  resulting  in  modified 
equations  (6)  and  (7)  which  include  an  effective  coupling  length,  instead  of  the  physical  length,  as 
suggested  in  [12].  The  effective  coupling  length  can  be  defined  as: 


Lt}} 

/\(f> 


L  + 


A(/> 


[/3^e{z)  -  (3zo{z)]dz 


(8) 

(9) 
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where  the  integration  is  i)erforniecl  along  the  axial  diroction  t)f  the  couph’r.  Cn  corresponds  to 
the  junction  between  the  coupled  guide  and  tlu-  connecting  arm.  while  is  chosc'ii  to  Ix'  some 
value  of  r  beyond  which  the  coupling  between  the  arms  is  negligible.  Due  to  th('  assutiied  ideal 
conditions  of  operation,  equations  (6)  and  (7)  cannot  account  for  any  loss  in  the  coupler.  In  [3S]. 
these  ideal  designs  have  been  applied  to  the  design  of  several  siib-millimeter-wave  couplers  and 
their  predicted  coupler  geometries  were  analyzed  by  the  FDTD  technique.  The  predicted  and 
numerically  derived  results  were  compared  in  an  effort  to  understand  the  effect  of  radiation  loss  on 
the  coupler  performance.  In  this  report,  an  example  of  a  mm-wave  image  line  coupler  is  investigated 
for  validation  purposes  and  a  sub-mm  wave  LRDW'  coupler  with  a  design  frequency  of  500  GHz  is 
presented. 

For  validation  purposes,  a  distributed  coupler  made  of  image  dielectric  line  is  investigated.  Fig¬ 
ure  12  shows  |5'i2|  and  jS'isI  of  the  image  guide  coupler  compared  to  the  experimental  results  from 
[39].  Considering  the  fact  that  the  experimental  results  were  obtained  through  an  implicit  mea¬ 
surement  with  the  assumption  that  the  coupler  is  lossless,  qualitative  agreement  between  theory 
and  experiment  is  satisfactory.  The  magnitude  of  Sn  is  found  to  be  less  than  -30  dB  in  the  whole 
frequency  range  which  agrees  with  the  experimental  results  in  [39].  Figure  13  shows  the  amount  of 
radiated  power  from  such  a  coupler  which  further  verifies  that  radiation  losses  can  not  be  neglected. 
It  is  worth  mentioning  that  the  FDTD  obtained  frequency  at  which  |<5'i2|  and  IS'iaj  intersect  agrees 
very  well  with  that  obtained  using  equations  (6)  and  (7)  in  conjunction  with  (8)  and  (9)  (see  Figure 
13  in  [39]).  It  should  be  mentioned  that  the  tapered  sections  of  the  dielectric  guide  are  modeled 
using  the  “staircase”  approximation. 

Using  (6)  and  (7)  in  conjunction  with  the  propagation  constants  of  the  even  and  odd  modes 
(obtained  using  2D-FDTD),  it  was  found  that  a  coupling  length  L  of  approximately  200^m  should 
give  a  3dB  coupler  with  a  design  frequency  of  500  GHz.  The  separation  between  the  two  LRDWs 
in  the  coupler  is  40  fim.  Figure  14  shows  the  scattering  parameters  of  this  coupler  with  a  coupling 
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Frequency  (GHz) 

Figure  12:  Magnitudes  of  S12  and  S13  of  the  image  guide  coupler  compared  to  the  experimental  results  from  [39].  S= 
2.4  mm,  image  guide  width=  4.8  mm,  image  guide  height=  2.4  mm,  Cr=  2.22,  L=59  mm,  6=26° ,  and  d=9.6  mm. 


Frequency  (GHz) 


Figure  13:  Radiation  loss  factor  (1  -  |SuP  —  15121^  —  |5'i3p  —  of  the  image  guide  coupler. 
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length  of  250/im.  The  sources  of  the  disrrepaiic}-  Ix'tweeii  th('S(‘  siimilaKHl  rc'sults  and  tlu'  ideal 
response,  which  should  have  a  :klB  design  fre(,iuenc\-  around  dOO  (ill/,  are  the  elfects  ol  t h<' junct ions 
and  the  coupling  between  the  connecting  guides.  .Morco\'er.  h  ig.  M  shows  radiation  loss  as  i)r('dicted 
bv  the  FDTD  method.  From  this  figure,  it  can  be  seen  that  at  490  Ciliz.  almost  of  the  input 


power  is  lost  as  parasitic  radiation  in  the  form  of  space  and  surface  waves. 

Figure  15  shows  the  results  for  the  same  coupler  with  a  coupling  length  L  of  300^/m.  In  this  case. 
Si2  and  Si3  intersect  at  488  GHz  with  a  value  of  -3.8  dB.  With  this  in  mind,  the  coupling  length  L 
was  increased  to  330//m  in  order  to  get  a  design  frequency  closer  to  500  GHz.  Figure  16  shows  the 
results  for  the  coupler.  In  this  case,  S12  and  S13  intersect  at  approximately  500  GHz  with  a  value 
of  -4  dB.  Finally,  a  coupler  with  GaAs  used  in  the  middle  layer  (i.e.,  £^=12. 85  instead  of  10)  was 
analyzed  and  its  response  was  found  to  be  similar  to  that  shown  in  Fig.  16  except  that  the  values 
of  S12  and  Si3  at  the  design  frequency  (500  GHz)  were  -5  dB  and  the  radiation  loss  was  somewhat 
higher  than  that  shown  in  Fig.  16. 


3.3  Matched  Termination 

Using  the  2D-FDTD  method,  it  was  found  that  in  the  LRDW  shown  in  Fig.  1,  most  of  the  power 
is  confined  in  the  substrate  just  underneath  the  ridge.  The  ridge  acts  only  as  a  guiding  element  for 
the  power  propagating  along  the  line.  ^Vith  this  fact  in  mind,  we  analyzed  an  open-end  LRDW  (by 
stopping  the  ridge  at  some  point  while  leaving  the  substrate  extended).  The  return  loss  of  such  an 
open-end  is  shown  in  Fig.  17.  A  return  loss  of  less  than  16  dB  in  the  frequency  range  450  to  550 
GHz  is  obtained.  Thus,  such  an  open-end  may  be  used  as  a  matching  termination  to  the  LRDW. 

Figure  18  shows  the  scattering  parameters  of  an  LRDW  coupler  with  port  4  open-ended.  It  can 
be  seen  that  S12  ^-nd  S13  intersect  at  505  GHz  with  value  of  -3.8  dB. 

Other  matching  mechanisms  can  be  used.  For  example,  a  lossy  material  that  absorbs  power 
at  sub-millimeter  waves  situated  at  the  end  of  the  LRDW  can  act  as  a  matched  termination. 
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Figure  14:  Scattering  parameters  and  Radiation  loss  factor 
:=40  fim,  L=250  ^m,  d=320  fim,  e=  26°. 
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of  the  LRDW  coupler  obtained  using  the  FDTD  technique 
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Figure  16:  Scattering  parameters  and  Radiation  loss  factor 
S=40  ^m,  L=330  fim,  d=320  /im,  6=  26^ 


Figure  17:  Return  loss  of  an  open-end  LRDW. 


Frequency  (GHz) 


Figure  18:  Scattering  parameters  of  the  LRDW  coupler  with  port  4  open-ended.  S=40  //m,  L=300  /rm, 


d=320 


0^  26". 
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Alternatively,  the  LRDW  may  be  transitioned  to  a  rectangular  waveguide  which  is  then  terminated 
properly. 

4  Design  Guidelines 

In  the  design  of  the  above  coupler,  the  following  steps  have  been  followed: 

1.  The  thicknesses  and  the  relative  permittivities  of  the  layers  are  specified. 

2.  The  2D-FDTD  method  is  used  to  obtain  the  dispersion  diagram  of  the  single  LRDW  under 
consideration.  This  is  done  in  order  to  know  the  value  of  the  effective  dielectric  constant  to  be 
used  at  the  front  and  back  absorbers  when  the  3D-FDTD  method  is  executed  to  analyze  the 
coupler. 

3.  The  separation  distance  between  the  two  coupled  LRDWs,  to  be  used  in  the  coupler,  is  speci¬ 
fied. 

4.  The  2D-FDTD  method  is  used  to  obtain  the  even  and  odd  propagation  constants  of  the  coupled 
LRDWs  structure  in  a  specific  frequency  range.  From  these  propagation  constants,  equations 
(6)  and  (7)  are  used  to  determine  the  coupling  length  needed  to  achieve  a  specific  coupling 
ratio  at  the  design  frequency. 

5.  The  3D-FDTD  method  is  used  to  analyze  the  coupler  with  the  above  specifications.  Most 
likely,  the  results  will  not  be  as  expected  because  of  the  parasitic  radiation  effects.  Thus,  the 
3D-FDTD  method  is  executed  for  different  coupling  lengths  until  the  desired  characteristics 
are  obtained. 
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5  Description  of  FDTD  Methods 


5.1  2D-FDTD 

The  2D-FDTD  method  is  based  on  the  technique  presented  in  [34.  40.  41].  In  tliis  technique,  tlie  3D 
standard  Yee’s  cell  [16]  is  reduced  to  a  2D  cell.  This  technique  can  be  used  to  obtain  the  dispersion 
diagram  of  any  waveguide  structure  or  transmission  line.  Thus,  instead  of  discretizing  a  3D  mesh, 
only  the  cross  section  of  the  line  under  consideration  is  discretized. 

Note  that  in  obtaining  the  dispersion  diagram  of  the  single  LRDW  shown  in  Fig.  1.  only  one 
sample  point  is  included  that  corresponds  to  /9=21000  rad/sec.  It  should  be  emphasized  that  in 
2D-FDTD,  the  propagation  constant  /?  is  chosen  first  and  then  the  frequencies,  at  which  there 
exist  modes  propagating  with  this  specific  /3,  are  obtained.  This  is  in  contrast  to  frequency  domain 
techniques  in  which  the  opposite  usually  happens.  The  space  steps,  Ax  and  Ay  are  carefully  chosen 
such  that  integral  numbers  of  them  can  approximate  the  various  dimensions  of  the  structure.  As  a 
rule  of  thumb  and  in  order  to  reduce  the  truncation  and  grid  dispersion  errors,  the  maximum  step 
size  is  chosen  to  be  less  than  1  /20  of  the  smallest  wavelength  occurring  at  the  highest  frequency  of 
interest.  Then,  the  Courant  stability  criterion  is  used  to  select  the  time  step  to  insure  numerical 
stability.  For  the  LRDW  under  study,  the  following  parameters  are  used:  Aa:=5.0  ym,  Ay=10  y.m 
and  A^=0.01489  ps.  The  number  of  time  steps  (ts)  is  chosen  large  enough  (e.g.,  10000)  to  let  the 
fields  reach  the  steady  state.  One  might  experiment  with  this  parameter  to  check  the  convergence 
of  the  final  result.  In  order  to  excite  the  computational  domain,  several  impulses  are  assumed  to 
exist  at  arbitrary  nodes.  This  will  give  all  modes  that  can  be  supported  by  the  guide.  However, 
one  may  change  the  program  such  that  a  specific  initial  condition  with  a  specific  kind  of  symmetry 
is  assumed.  It  is  also  necessary  to  specify  the  number  of  different  regions  that  the  computational 
domain  is  composed  of.  Then,  the  relative  permittivities  of  all  regions  have  to  be  entered  (notice 
that  anisotropic  materials  are  assumed).  Finally,  the  location  of  any  horizontal  or  vertical  perfectly 
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Figure  19:  Example  of  Fourier  Transform  Results 

electric  conductors  (PECs),  other  than  the  ground  plane,  have  to  be  specified.  Clearly,  there  are  no 
PECs  for  the  case  of  the  LRDW.  Once  the  FDTD  time  loop  is  done,  data  is  created  which  contains 
the  value  of  the  probed  field  component  as  a  function  of  time  step.  Then,  the  frequency  range  in 
which  the  discrete  Fourier  transform  of  the  probed  field  is  required  is  requested.  The  results  of 
this  Fourier  transform  contain  the  amplitude  as  a  function  of  frequency.  Fig.  19  shows  a  plot  of 
a  typical  example.  The  two  peaks  seen  in  the  figure  correspond  to  the  propagating  modes  that 
can  be  supported  by  the  LRDW  with  the  appriori  chosen  {3.  The  first  peak  corresponds  to  the 
dominant  mode  while  the  second  one  corresponds  to  the  first  higher  order  mode.  Finally,  choosing 
other  values  of  [3  and  obtaining  the  corresponding  frequencies  gives  the  dispersion  diagram  shown 
in  Fig.  8. 

Similarly,  one  can  obtain  the  dispersion  characteristics  of  the  even  and  odd  modes  of  the  coupled 
LRDWs.  However,  in  this  case,  one  needs  to  enforce  the  symmetry  to  be  able  to  tell  which  mode 
corresponds  to  the  odd  mode  and  which  one  corresponds  to  the  even  mode.  A  somewhat  modified 
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2D-FD7  D  molhod  is  used  in  th('  case  cjf  coupled  di('l('ct rics  to  ohtaiii  propagation.  1  lu'se  ju'oiuiga- 
tion  constants  will  be  used  later  in  conjunction  with  ('(|uations  (i  and  7  tcj  olitain  an  ('st iinat(' of  t he 
desired  coupling  length.  It  should  be  emphasized  that  in  tlu'  case  of  coupled  dieh'ct rics.  a  spc'cilic 
initial  field  distribution  is  used  to  enforce  the  ('ven  or  odd  symmetry. 

5.2  3D-FDTD 

The  3D-FDTD  method  is  based  on  the  conventional  Yee's  FDTD  scheme  [16].  It  is  necessary  to  run 
an  LRDW  through  line  before  analyzing  the  coupler  in  order  to  obtain  the  incident  pulse  waveform 
which  is  used  later  to  derive  the  scattering  parameters  as  was  discussed  in  the  theory  section.  It 
should  be  noted  that  the  definition  of  the  regions  here  is  similar  to  that  applicable  above  except 
for  the  fact  that  the  regions  are  three  dimensional  in  this  case  instead  of  being  two  dimensional. 
Figure  20  shows  the  field  component  at  some  distance  along  the  line  as  a  function  of  time  step. 
Some  reflections  can  be  noticed  due  to  the  imperfections  in  the  absorber  at  the  back  wall  of  the 
computational  domain.  This  is  due  to  the  fact  that  a  wide  band  pulse  (350-600  GHz)  was  used  as 
an  excitation.  Figure  21  shows  E^;  as  a  function  of  time  for  a  narrower  band  pulse  that  extends 
from  450-550  GHz.  In  this  case,  an  Ce//  of  4.0  is  used  at  the  front  and  back  walls.  It  can  be  noticed 
that  the  reflections  from  the  absorber  have  diminished  to  a  negligible  value. 

When  a  through  line  is  analyzed,  one  can  easily  get  the  effective  dielectric  constant  of  the  line 
as  a  function  of  frequency  using  transmission  line  theory  [18].  Two  probes  at  a  distance  of  lOAz 
are  usually  used  to  derive  the  propagation  constant. 

For  the  analysis  of  the  LRDW  coupler  (see  for  example  Fig.  18)  the  tapered  sections  have  to 
be  discretized  properly  using  the  “staircase”  approximation.  The  field  component  Ea,  is  probed  in 
4  locations,  each  at  a  different  port.  The  scattering  parameters  Sn,  S12,  813  and  S14  can  then  be 
calulated. 
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Pulse  Propagation 
LRDW 


Figure  21:  A  plot  of  the  field  components  Ex 
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OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Material 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence, 
reliability  science,  electro-magnetic  technology,  photonics,  signal 
processing,  and  computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


